Transport in chaotic quantum dots: 
effects of spatial symmetries which interchange the leads 
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We investigate the effect of spatial symmetries on phase coherent electronic transport through 
chaotic quantum dots. For systems which have a spatial symmetry that interchanges the source and 
drain leads, we find in the framework of random matrix theory that the density of the transmission 
eigenvalues is independent of the number of channels N in the leads. As a consequence, the weak 
localization correction to the conductance vanishes in these systems, and the shot noise suppression 
factor F is independent of TV. We confirm this prediction by means of numerical calculations for 
stadium billiards with various lead geometries. These calculations also uncover transport signatures 
of partially preserved symmetries. 
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I. INTRODUCTION 

Over the past years transport experiments on phase co- 
herent mesoscopic systems have attained unprecedented 
levels of sophisticationi 2 ^ It is now possible to precisely 
design and control the geometries of quantum dots while 
reducing the effects of impurity scattering to a level where 
the transport can be considered as purely ballistic (with 
scattering only off the confining boundaries) For 
geometries which give rise to chaotic classical motion 
universal system properties are expected^ 3 . Prototypi- 
cal chaotic cavities (such as the stadium and the Sinai 
billiard) do however feature geometrical symmetries that 
leave signatures on the transport properties j&LHiiL which 
we investigate in the present communication. In partic- 
ular, we consider systems with a lead-transposing reflec- 
tion symmetry which interchanges the source and drain 
leads that couple the dot to the electronic reservoirs [see, 
e.g., Figs. da) and (b)], and contrast them with systems 
that do not possess such a symmetry [see, e.g., Figs.^c) 
and (d)]. 

Our investigation is based on random-matrix theory 
(RMT) for the scattering matrixiiMi In fact, the present 
RMT problem is long-standing, dating back to almost 
10 years ago when first t heor etical and numerical results 
were presented in Refs. 1(1 7l respectively. These earlier 
works identified the correct invariant measure of the scat- 
tering matrix and they studied several statistical prop- 
erties of the conductance. Here we give the complete 
solution of this problem by deriving the joint probability 
density of transmission eigenvalues, which determine all 
stationary transport properties, for an arbitrary number 
of open channels supported by the leads attached to the 
quantum dot. The density of the transmission eigenval- 
ues turns out to take a particularly simple form. A strik- 
ing signature of the solution is the absence of any non 
trivial dependence of the ensemble-averaged conductance 




0.2 0.4 0.6 0.8 1 
T 



FIG. 1: Distribution of transmission eigenvalues p(r) (nor- 
malized to the sample size), N = 1, in stadium geometries 
with injection from the left: lead-transposing left-right sym- 
metry (a, b) and without global symmetry (c, d). The nu- 
merical data (histogram) is compared with the corresponding 
RMT- predictions (solid lines). Due to the horizontal align- 
ment of the entrance lead in (c) a fraction of back-reflected 
paths have an up-down reflected partner trajectory of equal 
length (see, e.g., pair of dashed lines). Although for a class 
of trajectories the corresponding partner is cut short by the 
exit lead (see, e.g., pair of solid lines) transport is strongly in- 
fluenced by this partial symmetry (see text). Cavity area (a) 
A = (4 + tt)/2, (b)-(d) A = 4 + tt, lead widths (a) d = 0.125, 
(b)-(d) d — 0.25, and entrance lead tilt angle in (d) S = 11.8°. 



and shot noise on the number of open transport channels. 
This prediction is confirmed numerically for stadium bil- 
liards with different geometries. Our theory also explains 
the deviations from standard random-matrix theory ob- 
served in earlier numerical investigations. 12 
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The random-matrix theory of transport 10,11 is based 
on Landauer's scattering approach, which describes the 
transport properties by the scattering matrix 
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composed of amplitudes for transmission (t, t') and reflec- 
tion (r, r') between the channels in the entrance and exit 
lead, respectively. We assume that both leads support 
the same number of open channels N. We investigate 
the transport properties in terms of the eigenvalues r„ of 
tH which determine the conductance Gn = G Yl 1 r n 
(where G — 2e 2 /h is the conductance quantum) and the 
shot noise power— Pn = (2eGV) J2 n=1 T n {X — r n ) (where 
V is the applied voltage) . For uncorrelated electrons the 
shot noise power is given by the Poisson value Pn — 2el, 
where I is the time averaged current. Fermi statistics, 
however, induces electronic correlations and deviations 
from the Poisson value which are customarily quantified 
by the Fano factor 



(P 



(P 



(En=l T «) 



(2) 



where (. . .) stands for ensemble or energy average. 



This gives an ensemble-averaged conductance (Gn) — 
Goo + 0(N°), where G m = NG/2, and a Fano factor 
F N = Foe + OiN- 1 ), where = 1/4. 

The bimodal distribution Eq. is only valid asymp- 
totically for large N. The leading finite-iV correction 
for the conductance is the well-known weak-localization 
correction (G N ) -Goo = -G/4 + O(iV~ 1 ) 1 i . Similar cor- 
rections exist also for other transport properties such as 
the Fano factor. Generally, they can be related to devia- 
tions of Pn{ t ) from Poo( T )j an d are most pronounced for 
a small number of channels. 

For a single channel (N — 1), the density of transmis- 
sion eigenvalues is given by^^ 
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hence the conductance and shot noise power is given by 
(Gi) /G = 1/3, (Pi) /(2eGV) = 2/15, respectively. The 
Fano factor is then F\ = 2/5. 

For two channels, from Eq. J3J, we find 



p 2 (r) =4t-3V7+^. 

\ T 
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Thus (G 2 ) /G = 4/5, (P 2 ) /(2eGV) = 9/35, and F 2 = 
9/28. These results show that for a small number of 
channels the corrections to the large- N asymptotics are 
clearly noticeable. 



II. JOINT PROBABILITY DISTRIBUTION OF 
TRANSMISSION EIGENVALUES 

Throughout this paper we assume the absence of mag- 
netic fields and spin-orbit scattering, hence, we consider 
time-reversal symmetric systems without Kramers de- 
generacy. 



A. Chaotic quantum dots without spatial 
symmetries 

As a benchmark for our results we use the well- 
established random-matrix theory for such systems in 
absence of any spatial symmetriesi°*ii (see the geome- 
tries in the bottom row of Fig.^l. The scattering matrix 
is then a member of Dyson's circular orthogonal ensem- 
ble (COE)jii and the joint probability distribution of the 
transmission eigenvalues r„ is given by 



W({r}) = [] 
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For a large number of channels N — > oo, the density of 
transmission eigenvalues pn(t) = (j2n ~ r ™)) a P" 
proaches the bimodal distributioni£*i& 



Pn(t) « Poo(t) = N 
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B. Chaotic quantum dots with spatial symmetries 

We now turn to systems with a spatial reflection sym- 
metry which interchanges the leads, such as the geome- 
tries of Fig.n](a) and (b). 

For a left-right symmetric structure, Fig. the scat- 
tering matrix has the structure 



S = 



r t 
t r 



(7) 



where r and t are symmetric TV x N matrices. Matrices 
of this structure can be cast into a block diagonal form 
by using the rotation matrix 6 



R = 
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S = RSR = 



si 
s 2 



(8) 



(9) 



where s\ = r + t and s 2 = r — t are symmetric unitary 
matrices. In terms of these matrices, the transmission 
matrix is given by t = ^(si — s 2 ) and the reflection matrix 
is given by r = \(s\ + s 2 ). 

The transmission eigenvalues are obtained from the 
matrix 



it 



2 — SlS 2 — (SlSrj) 
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which involves the unitary matrix Q = s\s\. Note that 
Q and can be diagonalized simultaneously. Hence, 
the eigenvalues exp(i0„) of Q determine the transmission 
eigenvalues by r„ = sin (0„/2). 

In random-matrix theory the matrices Si and S2 are 
assumed to be independent members of the COE. We 
then can show that the joint probability distribution of 
the eigenvalues exp(z6*„) is also given by the COE>i2ii^ 
The matrix Q can be symmetrized by the unitary trans- 
formation Q' = s 2 l ^ 2 Qsl/ 2 = s 2 X ^ 2 s\s 2 1 ^ 2 , which leaves 
the eigenvalues invariant. Moreover, the COE is invariant 
under the automorphism Q' — > U T Q'U, where U is an ar- 
bitrary unitary matrix which we identify with U — s 2 ■ 
Hence the eigenvalues of Q inherit the COE statistics of 
s%, and the joint distribution function of transmission 
eigenvalues takes the form 

w({r}) = n\[T k (i-T k )r^ ) 

k I CT fc =±lJ 
X Y[ | CT m\/ r m( 1 ~ T n) ~ ^n\J T n {\ - T m ) . 

(11) 

The probability density pn(t) can be obtained directly 
from the uniform distribution of cigcnphascs 9 n in the 
COE. This yields /9jv(t) = Poo(r) [given in Eq. (@J] ex- 
actly, for any value of N. Hence the ensemble-averaged 
conductance is given by (Gn) — Goo — NG/2, i.e., 
the weak-localization correction is absent. Moreover the 
Fano factor is given by Fn = Fqo = 1/4. The absence 
of any iV-dependent corrections is in striking contrast to 
the previously discussed case of asymmetric systems. 



C. Chaotic quantum dots with fourfold symmetries 

The results above also apply to systems with a 180° 
rotational symmetry mapping the leads onto each other. 
They can also be extended to incorporate further spatial 
symmetries which may be present in addition to the lead- 
transposing symmetry. For systems with two leads, the 
only remaining case is the fourfold symmetry as indicated 
in Fig. In this case, the S matrix can be written as 

s ={ S ol)> W 

where S e and S are the scattering matrices associated 
with the even and odd channels, respectively. Each ma- 
trix 5* e ( ) has the structure of Eq. (JJJ; thus each of them 
can be analyzed by the same procedure as in the previous 
lead-transposing symmetry case and the statistical trans- 
port properties are given by the superposition of the even 
and odd subsystems. For a cavity with a fourfold sym- 
metry we therefore have (Gn) = GN/2 and a constant 
Fano factor Fn = 1/4. 



III. NUMERICAL SIMULATIONS 

We now compare our theoretical results to numerical 
simulations of stadium billiards which feature the four 
different setups depicted in Fig. ^ The numerical data 
is obtained by solving the Schrodinger equation using a 
modular recursive Green's function method which allows 
for an efficient calculation of the scattering matrix even 
for a large number of open channels^ 

The most drastic effects of the lead-transposing re- 
flection symmetry are expected for the case of a single 
channel in the leads (N — 1). Figure ^ compares the 
numerical probability distribution p\ (t) with the analyt- 
ical predictions for the systems with and without a lead- 
transposing reflection symmetry (see upper and lower 
panel, respectively). Note that for each of the two pairs 
of geometries the numerical results (i) show a clear signa- 
ture of the absence/presence of the symmetry and (ii) we 
can see a good agreement with the analytical predictions. 

Several earlier studies have demonstrated that dynam- 
ical signatures of geometries can also be identified in 
terms of the Fano factori 12 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 Figure 
I2K shows the Fano factor for the symmetric geometries 
as a function of the number of channels. For the two 
systems with lead-transposing reflection symmetry our 
numerical simulations show an overall constant behavior 
of the Fano factor with the number of channels, around 
F = 1/4, in agreement with our modified random-matrix 
theory. The modes of the left-right symmetric cavity in 
Fig. 1(a) are identical to the modes with even index 
of the cavity in Fig. 1(b). Therefore, the presence of 
the odd channels in our fourfold symmetric cavity does 
not change the statistical properties of the transmission 
eigenvalues as predicted above. The flat behavior of the 
Fano factor in the symmetric cavities is nicely contrasted 
by the results for the two geometries without symme- 
try, Fig. 133: As predicted by the conventional random- 
matrix theory (absence of spatial symmetries) the Fano 
factor decreases as N increases and approaches its uni- 
versal value 1/4 for large 

Note, however, that in Fig. [2b the Fano factor of the ge- 
ometry with the horizontally attached lead (Fig.[TJ;) con- 
sistently lies above the RMT prediction for small channel 
numbers. This feature can be attributed to the fact that 
in this geometry a significant fraction of back-reflected 
classical tractories come in pairs of two equally long paths 
which are related to each other by an up-down mirror re- 
flection. This feature does, however, not apply to those 
trajectories the mirror images of which are cut short by 
the exit lead (see the two trajectory pairs in the inset 
of Fig.^J;). Since this "partial" symmetry enhances the 
constructive interference of reflected paths but spares all 
transmitted trajectories, the weak localization contribu- 
tion to the conductance and the shot noise power are 
increased. By tilting the entrance lead slightly (as in 
Fig-El) this effect should be destroyed, resulting in trans- 
port statistics which follow the RMT-prediction. Check- 
ing with our numerical data in Fig. [3 we note that this 
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FIG. 2: Energy-averaged Fano factor in different mode in- 
tervals. Upper panel: Numerical data for systems with a 
lead-transposing symmetry (see Fig.^i,b) compared with the 
corresponding mode-independent random-matrix prediction 
(F = 1/4, solid line). Lower panel: Analogous plot for sys- 
tems without global symmetry (see Fig.U,d). The deviations 
to the RMT-result (solid line, numerically obtained from Eq. 
© for N > 2) for the geometry depicted in Fig. 0: can be ex- 
plained in terms of the partial symmetry of short trajectories 
in the structure. 



assessment is indeed nicely corroborated. Given this ex- 
cellent agreement, one might wonder why recently pro- 
posed corrections to RMT due to "noiseless scattering 



channels"^2Si2Si2S do not apply in the present cases. Due 
to the small ratio of lead width d to cavity area A in 
our geometries (see caption Fig. ^) the criterion for the 
emergence of such fully transmitted or reflected channels, 
N > (fcf \^4)^ 2 r^ is, however, not fulfilled in the energy 
regimes studied above (for the highest energy considered 
here: N = 7 and {kp^fA) 1 ' 2 w 15). 



IV. SUMMARY 

We have studied analytically and numerically the ef- 
fects of spatial symmetries on electronic transport prop- 
erties of ballistic lateral quantum dots, modeled by a 
quantum chaotic cavity. Especially, we considered ge- 
ometries with a symmetry which maps the two leads onto 
each other. For such systems, random-matrix theory can 
be solved exactly for an arbitrary number of channels N 
in each lead. We predict that finite- TV corrections are ab- 
sent for all transport quantities, such as the conductance 
and the shot noise. This is confirmed by our numerical 
simulations. We further explored the effects of partial 
symmetries on transport which, as we showed, can yield 
significant corrections to the random-matrix predictions. 
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